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We describe radiative corrections to bremsstrahlung and their application to high energy collider physics, 
focusing on the applications to luminosity measurement, fermion pair production and radiative return. We review 
the status of one loop radiative corrections in BHLUMI and the /CA3MC, including cross checks with newer results 
developed independently for radiative return. We outline a YFS-exponentiated approach to the Drell-Yan process 
for LHC physics, including a discussion of the relevant radiative corrections. 



RADIATIVE CORRECTIONS 
BHABHA SCATTERING 



TO 



The BHLUMI Monte Carlo (MC) program[T] 
was developed as a precision tool for calculat- 
ing the small-angle Bhabha luminosity process 
at SLC and LEP, and with continued develop- 
ment, it will continue to be a valuable tool meet- 
ing the requirements of a next-generation linear 
e+e" colhder, such as the proposed ILC. Central 
to this program's success was an exact treatment 
of the phase space for n photon bremsstrahlung. 
A YFS-exponentiationj2] procedure allows all IR 
singularities to be canceled exactly between real 
and virtual emission processes to all orders. The 
leading soft photon effects are exponentiated, and 
IR-finite YFS residuals are then calculated ex- 
actly to the order required to reach the desired 
precision level. 

BHLUMI attained a total error budget of 
0.061% for LEPl parameters and 0.122% for 
LEP2 parameters for a typical calorimetric de- 
tector scenario. [3] To assure this precision level, 
it was necessary to calculate the most important 
unimplemented effect in BHLUMI4.04, which was 
the next to leading-log (NLL) contribution to 
the two-photon radiative corrections. The dou- 
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ble real photon and single real plus virtual pho- 
ton corrections to the small-angle Bhabha scat- 
tering process were calculated exactly in refs. 
|4I5| When added to the known two-loop vir- 
tual correction 6 , these results showed that the 
0{a-^L) corrections enter at the 0.027% level for 
LEPl parameters and 0.04% level for LEP2 pa- 
rameters. Here, L = \n{\t\/m1) is the "large log- 
arithm" entering into a leading log expansion. 

Implementing these exact 0{a'^) results in BH- 
LUMI would eliminate these contributions to the 
error budget. The only remaining unimplemented 
0{a'^) radiative corrections would then be up- 
down interference effects in which two virtual 
photons are exchanged between the e+ and e~ 
line, which are suppressed at small angles, and 
nominally enter at the level of 0.004% or less for 
angles below 9° -17181 These contributions, repre- 
sented by diagrams of the type shown in Fig. 1, 
could thus be safely neglected for SLC and LEP 
physics. 

For ILC physics, where the goal is to reach 
0.01% in the small-angle Bhabha luminosity pro- 
cess, it is desirable to carefully check the magni- 
tude of the up-down interference terms, and to 
implement them if they turn out to be signifi- 
cant. A key ingredient in the comparison, the 
five-point box integral appearing in the second 
and fourth diagram in Fig. 1, has recently been 
provided by the Looptools 2.2 package. 9^ A num- 
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Figure 1. One- loop up-down interference dia- 
grams for i-channel positron line emission. The 
internal photons may also be replaced by a Z bo- 
son. 



ber of other 0{a^) calculations which have ap- 
peared recentlv |lUlll) should also provide valu- 
able insight into effects which may need to be im- 
plemented in the small-angle Bhabha calculation 
to reach ILC precision specifications. 



RADIATIVE CORRECTIONS 
FERMION PAIR PRODUCTION 



TO 



Fermion pair production plays a critical role 
in extracting precision electroweak physics from 
e+e" colliders. This process is calculated by 
another YFS-exponentiated MC program, the 
/C/CMC.jni Again, the photonic radiative correc- 
tions play an essential role in calculating the YFS 
residuals through order . These have been cal- 
culated exactly, including finite-mass corrections, 
for initial state and final state radiation. |13| 

Using helicity spinor techniaues. |14ll5] a con- 
cise and stable representation for the 0{a^) ini- 
tial or final state radiation amplitude has been 
obtained, including finite-mass corrections. The 
matrix element for hard photon initial-state emis- 
sion with one virtual photon may be expressed as 

Mr^'' = ^Mf^-^'\fo + hh + hh), (1) 

where A^^^^'-^'' is the tree- level matrix element for 
single hard photon emission, fi are scalar form 



factors and li are spinor factors defined in ref. 

m 

The single hard photon cross section is 
of particular interest in radiative return 
applications |16ll7ll8| . where initial state radi- 
ation is used to reduce the effective beam energy, 
allowing a fixed energy machine to probe a range 
of energies. A MC program PHOKHARA was 
developed to calculate radiative return at <i> and 
B factories. |19I2(J| The same radiative corrections 
are relevant for a high-energy e^e~ collider inves- 
tigating physics around the Z peak, for example, 
ft is therefore useful to compare the radiative 
corrections obtained for both the /C/CMC and 
PHOKHARA in detail. Both calculations claim 
the same level of exactness, including the same 
diagrams as well as electron mass corrections 
relevant for coUinear bremsstrahlung. 

We have compared the virtual corrections to 
initial state hard-photon emission calculated in 
ref. (KR) for PHOKHARA to those cal- 

culated in ref. [T3| (JMWY) for the /C/CMC in 
the case of muon pair production. Analytically, it 
was found that in the absence of mass corrections, 
both expressions agree to NLL order (©(a^L) in 
the integrated cross section) . A compact ex- 
pression for NLL limit of the matrix element was 
obtained in ref. I13J , where it was shown that the 
terms /i and /2 in eqn. vanish to NLL order, 
and the helicity-averaged NLL limit of /o is 



Re (/r^) 



27rReBYFs(s) + L-l 
21nri ln(f - rz) - ln^(l - ri) 
ri(l - ri) 



+ 31n(f-ri) 
+ 2Sp(ri) + (ri 



1 



(2) 



where L = hi{s/m'^) is the "large logarithm" in 
the leading log expansion, ri = 2pi -k/s measures 
the inner product of one of the incoming fermion 
momenta pi with the hard photon momentum fc, 
Sp(a;) — Li2(x) is the dilogarithm (Spence) func- 
tion, and 
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is the IR-divergent virtual YFS form factor. 

Since the two results are known to agree with 
the NLL limit calculated using eqn. , the NLL 
limit is subtracted in each case, permitting the 
NNLL contributions and collinear mass correc- 
tions to be investigated in the context of the 
/C/CMC. Fig. 2 shows the resuhs of a /C/CMC 
run calculating the NNLL contribution to muon 
pair production at a CMS energy of 500 GeV for 
10^ events, both with and without the mass cor- 
rections. The cross-section is integrated up to 
a radiated photon energy fraction of Umax (with 
V — ri + r2) using the YFS residual (3l for 
one hard photon at 0(0^), subtracting the NLL 
contribution obtained using eqn. The result 
is normalized with respect to the non-radiative 
Born cross section for muon pair production. 
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Figure 2. Comparison of NNLL contributions to 
the virtual correction to the single hard photon 
cross section for muon pair production at a CMS 
energy of 500 GeV. 



It is found that the maximum difference be- 
tween the complete KR and JMWY results (from 
the next-to-last bin) is 1.6 x lO"^ units of the 
Born cross section. Most of this apparently 
comes from differences in the treatment of the 



mass corrections. KR uses an expansion in 
m\/riS, while JMWY uses a technique developed 
by Berends et a/. [21] for adding the mass cor- 
rections required in collinear limits to a calcu- 
lation obtained using massless spinors. With- 
out mass corrections (comparing the massless 
points), the results agree to within a part per 
million. This agreement is better than noted 
previouslv |23l25l26l27| due to improvements in 
the stability of the algorithms used. Direct com- 
parisons of the PHOKHARA and KKMC pro- 
grams have also been conducted. |28I29| 

It is interesting to note that the size of the 
NNLL part of the corrections implemented by 
JMWY never exceed 4 x 10"^ up to u^ax = 0.975, 
and reach —3.25 x 10"^ in the last bin, where 
■f^max = 0.9875. This suggests that for most pur- 
poses, the considerably simpler NLL result repre- 
sented by eqn. (jSJ will suffice. 

3. THE DRELL-YAN PROCESS 

The Drell-Yan process plays a role at hadron 
colliders which is as basic as the Bhabha scatter- 
ing or pair production cross section at e^e~ col- 
liders. In fact, W and Z production has been pro- 
posed as the luminosity process for the LHCpU] 
A fully exclusive calculation of the parton-level 
cross sections is needed at the 1 — 2% level for 
upcoming LHC physics. These cross sections are 
currently known at the 10% level, using NLO ma- 
trix elements. 

While NNLO results are available for 
the integrated cross section "ST and rapidity 
distribution 32 , a fully-exclusive NNLO cross 
section needed for a MC event generator is not 
yet available. Moreover, electroweak radiative 
corrections will be required as well. Reaching 
the desired LHC precision will require correc- 
tions of order a1 and order aew, including mixed 
©(ofsCkew) contributions. Examples of the latter 
diagrams are shown in Fig. 3. 

The hard parton-level processes must be cal- 
culated and combined with PDFs in a MC pro- 
gram designed to generate the desired distri- 
bution of partons plus mixed QCD and QED 
bremsstrahlung. This will require a careful im- 
plementation of the multiple gluon and photon 
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Figure 3. Examples of mixed photonic and glu- 
onic initial-state radiative corrections: (a) real 
photon + gluon emission, (b) real photon with 
QCD loop, (c) real gluon with electroweak loop, 
(d) 2 loops, mixed QCD/electroweak. The final 
state fermions from the decay are not shown. 



phase space. Experience in the electroweak sec- 
tor suggests that YFS exponentiation will pro- 
vide a strong tool for implementing the multiple- 
emission phase space, giving very precise control 
over the soft and coUinear limits. 

The large numbers of diagrams creates a chal- 
lenge for obtaining an expression that can be 
evaluated quickly enough for MC implementa- 
tion, and evaluated in a numerically stable man- 
ner. Common reduction methods based on 
the Passarino-Veltman technique Plj can produce 
millions of terms, which are both slow to evaluate 
in a MC setting, and prone to numerical insta- 
bilities due to the large numbers of terms added 
and potential cancellations among them. Thus, 
a significant part of this problem will involve de- 
veloping and testing new methods for organizing 
and calculating the terms in a stable manner. 

Once the parton-level matrix element is ob- 
tained, it may be incorporated into a "QCED- 
Exponentiated" MC program, implementing a 
procedure similar to YFS exponentiation in a 
combined QED and QCD setting to construct 
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the exact phase space for multiple gluon and 
photon radiation. This requires extending the 
YFS calculus to non-abelian gauge theory, with 
due care in handling the genuine non-abelian sin- 
gularities which arise in QCD. The QCD and 
QED exponentiation will be conducted at order 
a^L on an event-by-event basis in the presence 
of showers without double-counting of soft and 
coUinear emission effects. Further details of this 
construction will be presented elsewhere in these 
proceedings. j33j 



4. SUMMARY AND OUTLOOK 

A careful calculation of higher-order brems- 
strahlung corrections led to a precision tool (BH- 
LUMI) for Bhabha luminosity calculations. In- 
corporating the second order photonic corrections 
obtained to test BHLUMFs precision may be 
enough to reach the 0.01% level proposed for the 
ILC. Calculating the O(a^) up -down interference 
contribution will help to clarify this. A number 
of second-order Bhabha scattering results are now 
available which should be useful for testing BH- 
LUMI's precision at this higher level. 

We have also described cross-checks on the 
second-order photonic radiative corrections for 
fermion pair production developed for the 
/C/CMC. Comparisons with similar initial-state ra- 
diative corrections developed for PHOKHARA 
show agreement on the order of 10~^ in units of 
the Born cross section at ILC energies. Remain- 
ing differences may be attributed primarily to dif- 
ferences in the handling of finite mass corrections 
in the coUinear limits, and possible residual nu- 
merical instabilities for photons radiated near the 
fermion pair production threshold. 

Finally, we have outlined a program for carry- 
ing the successes of the YFS-exponentiatcd MC 
framework developed for SLC and LEP physics 
into the hadronic realm of the LHC, where a pre- 
cise calculation of the parton-level diagrams can 
be combined with a combined QCD and QED ex- 
ponentiation framework, with complete control of 
the multiple photon and gluon phase space, to de- 
velop a MC for the DrcU-Yan process which can 
reach the LHC precision requirements. 
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